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Abstract 

Formation control is concerned with the design of control laws that stabilize agents at given distances 
from each other, with the constraint that an agent's dynamics can depend only on a subset of other 
agents. When the information flow graph of the system, which encodes this dependency, is acyclic, 
simple control laws are known to globally stabilize the system, save for a set of measure zero of initial 
conditions. The situation has proven to be more complex when the graph contains cycles; in fact, with 
the exception of the cyclic formation with three agents, which is stabilized with laws similar to the ones 
of the acyclic case, very little is known about formations with cycles. Moreover, all of the control laws 
used in the acyclic case fail at stabilizing more complex cyclic formations. In this paper, we explain 
why this is the case and show that a large class of planar formations with cycles cannot be globally 
stabilized, even up to sets of measure zero of initial conditions. The approach rests on relating the 
information flow to singularities in the dynamics of formations. These singularities are in turn shown 
to make the existence of stable configurations that do not satisfy the prescribed edge lengths generic. 

Index Terms 

Formation control, Decentralized control, Global stability, Bifurcations, Singularities. 

I. INTRODUCTION 



Consider the following problem, depicted in Figure la Three autonomous agents with posi- 
tions Xi,x 2 and x 3 evolve in the plane according to first order dynamics, agent 1 observes the 
position of agent 2, agent 2 the position of agent 3 and agent 3 the position of agent I. Can the 
agents stabilize at prescribed distances di,d 2 ,d 3 from each other? 
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Problems of this type, which fall under the broader class of decentralized control problems, 
have been a focal point of attention of control theory for the past decade or more, as they arise 
in a wide variety of natural (think schooling, herding, etc.) and engineering situations (think 
autonomous vehicles or decentralized power systems). The vast majority of control laws proposed 
in these contexts are so-called gradient control laws, named after the fact that every agent tries to 
minimize its own objective function. While such control laws work well in centralized systems 
with a small number of agents, they become ineffectual at stabilizing a given configuration 
when this number increases or when the system is decentralized due to the appearance of a 
very large number of stable configurations that do not respect the desired inter-agent distances. 
The objective of this paper is to address the effect of decentralization on the appearance of 
stable, undesired equilibrium configurations. In order to do so, we will consider a broader class 
of control laws than gradient laws and analyze a particular four-agents formation, called the 
2-cycles. This four- agents formation was exhibited in [OQ to illustrate the shortcomings of the 
current methods in formation control and decentralized systems, as it resisted attempts to either 
define globally stabilizing control laws or prove their non-existence. 

We will use ideas from singularity and bifurcation theory O, O to show that the 2-cycles, 
and trivially systems containing it as a subformation, are not globally stabilizable. Singularities 
have not often appeared in the study of global control design, the reason behind this fact is that 
they are, in general, easily avoided by considering a small perturbation of the system and can 
thus be made irrelevant to the dynamics. By opposition, we will show here that the information 
flow constraints inherent to decentralized control can make such singularities unavoidable. 



Fig. 1: (a). Three agent in a cyclic formation in the plane. Agent 1 observes agent 2, which observes agent 3 
which in turn observes agent 1. An arrow pointing from agent i to agent j thus implies that the dynamics of agent 
i is allowed to depend to the state of agent j. (b). The 2-cycles formation. 

Going back to the example of three agents in the plane, there are at least two ways in which 
such systems arise in practical situations: 




;r 2 
(a) 



(b) 
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• Every agent chooses independently the distance di at which to stabilize from another agent. 
Can we design control laws so that when the agents evolve, they will all stabilize at their 
desired distances from almost all initial positions? 

• A central authority decides on a target configuration and sends each agent only part of the 
complete description of the target, e.g. by sending them only the inter-agent distance they 
have to satisfy. Can we design a control law that will stabilize every target configuration? 

These two points of view of course lead to the same mathematical problem, though they arise 
in different contexts. In the first case, the question may arise in the study of natural flocks or non- 
cooperative settings, where agents are incapable or not wanting to communicate their objectives. 
In the second case, scenarios involving UAV's which are designed to achieve a global objective, 
but which for secrecy or efficiency reasons are only given part of the global objective, can 
give rise to such questions. For a more detailed discussion of local and global objectives of a 
decentralized system, see [|4]|. 

We now summarize the extant relevant work and methods used, see [51 for more details. We 
first mention that most issues in formation control arise when the information flow graph of 
the system contains cycles. Indeed, a cycle-free formation, which necessarily has a leader (i.e. 
an agent that does not observe any other agents) can easily be handled via the use of different 
time-scales depending on how distant from a leader an agent is. See jH for a detailed sketch 
of how a gradient-based control law globally stabilizes, save for a set of initial conditions of 
measure zero, acyclic formations. For related wor, dealing with undirected formations, we refer 
to and references therein. 

Remark 1. We mention here that the existence of a set of initial conditions that will not lead 
to the desired configuration stems from the nontrivial topology of the state-space of formation 
control, essentially the topology of a complex projective space [8]. We say that a control law 
almost surely stabilizes an equilibrium if the system stabilizes that equilibrium from almost all 
initial conditions. We come back to this in Section 

Local stabilization of formations that contain cycles has similarly been investigated in [|9]|, 
where it is shown that the widely-used gradient law can be modified by adjusting some gains 
to stabilize any desired target configuration. One caveat to the result is that the gains are 
not evaluated locally by the agents, resulting in a centralized "design phase" followed by a 
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decentralized "implementation phase" for the agents. It is shown in ifTOl that the gains cannot 
be locally evaluated by the agents in the case of the 2-cycles formation, i.e. that a decentralized 
design phase is not possible. 

This leaves us with the case of global stabilization of formations that contains cycles. The 
extant work consists of the thorough analysis of the triangle formation, or 1 -cycle, done in ifTTIl 
and related publications, where it is shown that the gradient control law almost surely globally 
stabilizes almost any configurations. We show in this paper that the second simplest formation 
with cycles, i.e. the 2-cycles formation, cannot be globally stabilized — even modulo sets of 
measure zero of initial conditions — by a broad class of control laws, including the control laws 
used in prior work. 

In order to prove the main result, we introduce two definitions. First, the already mentioned 
almost sure stability, in Section [n} which is needed to formalize the idea of global stability 
modulo sets of measure zero which has appeared in prior work on formation control. The 
second is the one of robustness for nonlinear systems; in a few words, we say that a control 
law is robust if its effect (e.g. stabilization) persists after small perturbations of the dynamics. 
This is closely related to structural stability of systems and requires some use of transversality 
conditions, discussed in the Appendix, to make rigorous. 

We conclude this introduction by an example: we illustrate how the gradient-based control 
law used in many works on directed formation control 0, [|9]|, IfTTTl fails to globally stabilize 
the 2-cycles formation by stabilizing around an undesired configuration. 

Let Xi E M 2 be the positions of the agents and di be positive numbers, the target edge lengths. 
The decentralized system is explicitly given by: 



ii = 


(11*1 


- x 2 \ 


- d 1 )(x 1 


- x 2 ) + (\\xi - x A \\ - d 5 )(x 1 - x A 


X 2 = 


(11*2 


~ 2*| 


- d 2 )(x 2 


~ 2*) 


x 3 = 


(11*3 


— X\\ 


- 4)(*3 


- x x ) 


24 = 


(11*4 


~ 2*| 


— C/ 4 )(X4 


- 2*) 



(1) 



V 



We show in Figure [2]the results of simulations for the vector of target distances d\ = 2.0, d 2 = 
2.6, d 3 = 2.0, g?4 = 1.4, d 5 = 3.3 and illustrate the appearance of a stable undesired configuration, 
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(a) Configuration Dl (b) Configuration D2 (c) Configuration Al 

Fig. 2: Simulation results for the decentralized system of Equation ([T]) with a d E C c . The dotted 
lines represent the trajectories followed by the agents. Both configurations D\ and D 2 represent 
design target frameworks. D 1 is locally stable whereas D 2 is locally unstable. The configuration 
Ai is also an equilibrium of the dynamics, though one for which the target distances are not 
respected: we call it ancillary equilibrium. A linearization of the system gives that the spectra 
of the Jacobians are given by (-17.5 + 1.3i, -17.5 - 1.3z, -11.9, -7.9, -0.6), (0.6, -18.6 + 
3i, -18.6 - 3i, -9.4 + 3.1z, -9.4 - 3.H) and (-23.4 + 4.8z, -23.4 - 4.8a, -11 + 2.8i, -11 - 
2.8i, —1.6) for the configurations Dl, D2, and Al respectively. Hence Al is locally exponentially 
stable and thus there is an open set of initial conditions that lead to an ancillary configuration. 

which is accompanied by an unstable desired configuration. 

II. Almost sure stability 

Consider the control system 

x = f(x,u(x)) (2) 

where x £ M, a smooth manifold, and all functions are assumed smooth. To justify the definition 
of almost sure stability, we restrict ourselves to the case of hyperbolic dynamics, i.e. having the 
property that the linearization of the dynamics at a zero has eigenvalues with non-zero real-parts. 

We are interested in global results about stabilization of an equilibrium configuration. From [[HI, 
we know that because of their invariance to rigid transformations of the plane, formation control 
systems evolve on the manifold M = CP(n — 2) x R, which is of dimension 2n — 3. The 
well-known Poincare-Hopf equality, which we illustrate in the case of the circle in Figure [3j 
relates the index of isolated zeros of differentiable vector fields on M to a global topological 
characteristic of M, its Euler characteristic lfi~2ll . 

On the one-hand, the Euler characteristic of M is known to be n — 1. For x an isolated zero 
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Fig. 3: We represent three continuous vector fields on the circle. The one on the left has no 
equilibrium, the one in the center has one stable equilibrium and one unstable equilibrium and 
the one on the right has two stable equilibria and two unstable equilibria; the indices of the 
equilibria are indicated on the figure. Because the sum of the indices is constrained to be zero, 
there is an even number of equilibria and, in particular, no continuous system on the circle can 
be globally stable. 



tells us that a stable equilibrium has an index of —1 (since the dimension of M is odd). Putting 
these two simple observations together, we conclude that global stabilization of an equilibrium 
is not possible for formation control, since stabilizing a zero will force the appearance of other 
zeroes ai to satisfy the Poincare-Hopf equality: 



On the other hand, we know from the Hartman-Grobman theorem that the dimensions of the 
stable and unstable manifolds of an equilibrium are given by the numbers of eigenvalues of the 
Jacobian with negative and positive real parts respectively. The attraction basin of an equilibrium 
is thus of codimension at least one (in other words, it is a thin set) unless all its eigenvalues 
have negative real parts. 

Hence, from a practical standpoint, if we only require that the control u(x) makes one equi- 
librium stable and all other equilibria either saddles or unstable, we obtain a system that behaves 
essentially like a globally stable system — since a vanishingly small perturbation would ensure that 
the system, if at a saddle or unstable equilibrium, evolves to the unique stable equilibrium — while 
leaving more room to possibly satisfy the Poincare-Hopf equality. We formalize and elaborate 



of a vector field, the useful formula 



ind(xo) = sign I det( 




ind(a;o) + 2, ind(aj) = n — 1 =>- >, ind(aj) 
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(a) (b) 

Fig. 4: Consider a control system defined on the sphere M. Assume that we want to design a 
control law such that the system stabilizes almost surely to either D\ or D 2 as depicted in (a) 
above. We call D\ and D 2 design equilibria and write Ed = {Di,D 2 }. A continuous feedback 
control law u(x) on the sphere that has either D x or D 2 , or both, as zeros may introduce additional 
zeros, which we call ancillary equilibria. We assume that there are two such equilibria A\ and 
A 2 and write E a = {Ai, A 2 }. The control law u(x) makes the system almost surely stable if at 
least one element of Ed is stable and no element in E a is stable. 



on this observation here. 

Let Ed be a finite subset of M containing configurations that we would like to stabilize via 
feedback. All configurations in Ed are equally appropriate for the stabilization purpose. We are 
thus interested in the design of a smooth feedback control u(x) that will stabilize the system to 
any point x G E&. We call these points the design targets or design equilibria: 

Ed = {xq G M s.t. x is a design equilibrium} 

Let 

E = {x G M s.t. f(x ,u(x )) = 0}, 

the set of equilibria of (|2]). We assume that £ is finite. 

As explained above, when the system evolves on a non-trivial manifold, the Poincare-Hopf 
equality, or the more refined Morse inequalities lfT3l . make it unreasonable to expect that there 
exists a control u(x) that makes the design equilibria the only equilibria of the system, i.e. such 
that Ed = £■ We call the additional equilibria, that are introduced by the non-trivial topology of 
the space, ancillary equilibria: 

E a = E — Ed- 
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We decompose the set £ into stable equilibria, by which we mean equilibria such that all the 
eigenvalues of the linearized system have a negative real part, and unstable equilibria, where 
at least one eigenvalue of the linearization has a positive real part. Observe that under this 
definition, saddle points are considered unstable. In summary: 

s S s u S u 

where 

£ s = {xo £ £ \ x is stable} and E u = {x 6 £ | x is unstable}. 

With these notions in mind, we introduce the following definition: 

Definition 1. Consider the smooth control system x = f(x,u(x)) where x G M and the set £ 
of equilibria of the system is finite. Let £d C M be a finite set. We say that £d is 

1) feasible if we can choose a smooth u(x) such that £d H £ ^ 0. 

2) almost-surely stabilizable if we can choose a smooth u(x) such that £ s C £d- 

When the set £d is clear from the context, we say that the system is feasible or almost-surely 
stable^] 

The set £ d is feasible if we can choose u(x) such that at least one equilibrium of the system 
is a design target. It is said to be almost surely stable if the system stabilizes to £d almost surely 
for all initial conditions on M. The usual notion of global stability is a particular instance of 
almost-sure stability; indeed, it corresponds to having u(x) such that £ d = £ = £ s . 

Looking at the contrapositive of this definition, a system is not almost-surely stable if there 
exists a set of initial conditions of codimension zero that leads to an ancillary equilibrium. 

Example 1. Consider a system 

x = x{l — kx 2 ) 

where k 6R is a feedback parameter to be chosen by the user. We show that any £ d C (0, oo) 
is not almost-surely stable. We first observe that the system has an equilibrium at and two 
equilibria at x — ± \/l/k if k > 0. The system is thus feasible for any £d C M. The Jacobian 

'in the earlier publication |14|, we termed almost-sure stability type- A stability 
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of the system is 1 at x = and —2 at x = ± ^/ljk. For k > 0, ?/ze above says that 

£ = {0, ±/lA} = { y^Tfc} U {0, -y/l/k} ■ 

From the linearization of the system, we have that 

S s = {±y^T/k} and S u = {0}. 

We conclude that £ s Sd and the system is not almost-surely stable. Indeed, all initial conditions 
Xq < result in the system stabilizing at an ancillary equilibrium. 

III. Genericity and robustness 

We now move on to the first of the two main technical ingredients necessary for the proof of 
the results below. The second ingredient, singularities and bifurcations, is presented in the next 
section. 

Let V be a binary-valued function on a topological space S, indicating whether a given 
property is satisfied. In more detail, if u G S, we say that u satisfies V is V(u) = 1. We have 
the following definition: 

Definition 2 (Robustness). An element u of a topological space S satisfies the property V 
robustly if for all u' in a neighborhood of u in S we have V{u') — 1. A property V is robust if 
there exists an open set U C S such that V{U) = 1. 

The property we will be dealing with here is stability: we want to find a u that stabilizes a 
system around an equilibrium, and desire the stabilization to be robust. In practical terms, if a 
property satisfied only at non-robust it's, then it fails to be satisfied under the slightest error in 
modelling or measurement. 

Related to robustness is the notion of genericity: a property V is generic for a topological 
space S if it is true on an everywhere dense intersection of open sets of S. 

Everywhere dense intersections of open sets are sometimes called residual sets fl2). 

Remark 2. We emphasize that when we seek a robust control law u(x) for stabilization, we 
seek a control law such that the equilibrium that is to be stabilized remains stable under small 
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perturbations in u{x). The equilibrium, however, may move in the state space. For example, 
assume that the system 

x = f(x, u(x)) 

has the origin as a stable equilibrium. If for all u(x) in an appropriate set of perturbations, the 
system 

x = f(x, u(x) + eu(x)) 

has a stable equilibrium at a point z(e) near the origin, then the control law u(x) is robust. If, 
on the contrary, the equilibrium disappears or becomes unstable, then u(x) is not robust. 

If 1 — V, the negation of V, is generic, then there is no robust u that satisfies V. Indeed, 
if 1 — V is generic, then V is verified on at most a nowhere dense closed set. In particular, 
V is not verified on an open set. The main tool to handle genericity are jet spaces and Thorn 
transversality theorem. We will use the results in some parts below and refer the reader to the 
appendix for more information. 

IV. Singularities, transcritical bifurcation and the logistic equation 

We recall a few definitions from dynamical systems theory. Consider a dynamical system of 
the form 

x = fn(x) (3) 

where x G M, an n— dimensional manifold, and p E M fc is a vector of parameters on which the 
system smoothly depends. 

Definition 3 (Hyperbolic and singular equilibria and bifurcation value). 

1 ) An equilibrium x is called hyperbolic if the eigenvalues of the linearization at x have 
non-zero real-parts. It is called singular or degenerate otherwise. 

2) A value po in the parameter space M. k for which the flow of ^ has a singular equilibrium 
is called a bifurcation value. 

A. The logisitc equation 

The logistic equation, which is often used to describe systems in which two competing 
effects — such as supply and demand or predator and prey — are at play, is the one-dimensional 
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Equilibria 




unstable 



stable 







Fig. 5: The logistic equation undergoes a transcritical bifurcation when /i = 0. The equilibrium x — is stable 
for jj, < and unstable for fx > 0. 



ODE given by 



This equation displays what is called a transcritical or transfer of stability bifurcation at ll = 0, 
which we explain here. Observe that the system has two equilibria, one at x = and one at 
x = jx, which coalesce when p, = 0. The linearization of the system about x is 



From this linearization, we see that for ll > 0, the equilibrium x = is unstable whereas 
the equilibrium x — p is stable. The situation is reversed for p < 0. We conclude that at the 
bifurcation value p = 0, the two equilibria coalesce and exchange their stability properties. We 



We show below that the 2-cycles behaves similarly to the logistic equation in the sense that 
they both exhibit the same type of singularities or bifurcation. The most common approach used 
to gain some understanding about the behavior of a dynamical system near a singularity relies 
on the use of the center manifold theorem |[T2l . This theorem establishes the existence of a 
nonlinear change of coordinates, valid near the equilibrium, where the dynamics can be put in 
a so-called normal form which is more amenable to analysis. The logistic equation as given 
in Equation [4] is such a normal form. This approach is without much hope for success for 
our purpose unless the control law u is fixed. Indeed, the change of variables involved in the 
analysis depends on the control u, and tracking the effect of this dependence through the whole 



x = x(/i — x). 



(4) 




{(i — x) — x — n — 2x. 



depict the above in Figure [5] This figure is to be compared to Figure 10 
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procedure is not feasible for broad classes of control laws. 

In order to sidestep this difficulty, we have recourse to the following result of Sotomayor [fT5l . 
which characterizes the generic behavior of dynamical systems near non-hyperbolic fixed-points 
without recourse to the center manifold. 

First, recall that for / : M 2 — > M 2 a twice differentiate function, its Jacobian is defined as 



df 
dx 



dfi dfi 

dx dy 

dh dh 

dx dy 



Assuming that the Jacobian has a zero eigenvalue, we denote by v and w corresponding right 
and left eigenvectors. The Hessian of/isa2x2x2 tensor with entries 

df\ dfj 
dx 2 J ijk dxjdx k 

Hence 



T 

W 



( v v \ _ ST w Ji v v 
dx 2 ' ^ % dxjdxk 3 

ijk J 



We have the following theorem: 



Theorem 1 (Sotomayor). Let x = f^x) be an ODE in W 1 depending on a scalar parameter \x, 
with f twice differentiable in x and \i. For p, = p , assume that the system has an equilibrium 
xo satisfying the following conditions: 

1) -j£r\ Xo has a unique zero eigenvalue with left and right eigenvectors w and v respectively. 
The other eigenvalues are negative. 

2) w T ^\ w v = 

3) w T %^\ X0 (v } v) and w T Qfc\ XOtfH) v ? 

Then the phase portrait is topologically equivalent to the phase portrait of the logistic equation, 
i. e. we have a transcritical bifurcation about Xq for p = fi . Thus around fi = /% there are two 
arcs of equilibria whose stability properties are exchanged when passing through p . Moreover, 
the set of equations x = fn(x) which satisfy conditions (1), (2) and (3) above is generic in the 
space of smooth one-parameter families of vector fields with an equilibrium at x , p with a 
zero eigenvalue. 
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V. Formation Control 

Let G = (V, E) be a graph with n vertices — that is V = {i>i, i>2, • • • , v n } is an ordered set of 
vertices and E C V x V is a set of edges. The graph is said to be directed if (i>j, i> j) € i£ does 
not imply that (vj,Vi) E -E. We let \E\ = m be the cardinality of We call the outvalence of 
a vertex the number of edges originating from this vertex. 

Directed graphs are used to encode the information flow in decentralized control problems. We 
follow the convention that an arrow leaving vertex Vi for vertex Vj means that agent % measures 
the relative position — relative to its own location — of agent j. 

Assume that the edges are ordered. The mixed-adjacency matrix of a graph G = (V, E) is a 
\E\ x |V| matrix whose entry (i, j) is —1 if edge e,- L originates from vertex vj, 1 if edge e, ends 
at vertex Vj and otherwise: 

Definition 4 (Mixed adjacency matrix). Given a directed graph G = (V, E), its mixed adjacency 
xm is defined by 



matrix A m E 



A — / 



-1 if ei = (vj,v s ),v s E V 
+1 ifei = (y k ,Vj),v k eV 







otherwise. 



The edge-adjacency matrix is a \E\ x \E\ matrix whose entry is —1 if edge e,- L and edge 
ej originate from the same vertex , 1 if edge ends at the vertex where edge ej starts and 
otherwise. Notice that A C:i j is zero if edge starts where edge ej ends and that the diagonal 
entries are —1: 



Definition 5 (Edge- adjacency matrix). Given a directed graph G = (V, E), its edge-adjacency 
matrix A e E R mxn is defined by 

-1 if e { = (v s , vt), ej = (v s , v f/ ,v s , v t , v t , E V 
4-e,ij = { 1 ifei = (v s ,vt),ej = (v t ,v s i),v s ,v s >,v t E V 
otherwise. 



We will often encounter the matrix A m ® / where ® is the Kronecker product and / the two- 
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Fig. 6: The 2-cycles formation. 



(2) 

by-two identity matrix. In order to keep the notation simple, we write A m for this Kronecker 
product. 

Example 2. The mixed-adjacency and edge-adjacency matrices of the 2-cycles of Figure [6] are 



A, 



-110 
0-110 
10-10 
1-1 

-10 1 



and A P 



-1 1 

-1 

1 


-1 



respectively, where edge i is labelled by Z{ as in Figure [6| 




1 

-1 
1 




-1 



1 
-1 

1 -1 



(5) 



□ 



A. Rigidity 

We briefly cover the fundamentals of rigidity and establish the relevant notation. We refer the 
reader to lTT6l for a more detailed presentation. We call a framework an embedding of a graph in 
M 2 endowed with the usual Euclidean distance, i.e. given G = (V,E), a framework p attached 
to a graph G is a mapping 

p : V R 2 . 

We write x, for p(vA. We define the distance function 5 of a framework with n vertices as 



8{p) : R 



. TO>2n ^ R n(n l)/2 . ^ _ _ ^ _ Jy^ _ _ _ _ ? _ ^ ^ y^ _ ^||2j ^ 



where IR + = [0, oo); i.e. 5(p) evaluates all the pairwise distances between edges. We denote by 
5(p)\e the restriction of the range of 5 to edges in E. 
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For a graph G with m edges, we define the set of feasible edge-lengths 

C = {d = (d u . . . , d m ) G for which Bp with 5(p(V))\ E = Vd} , 

where the square root of d is taken entry-wise. Properties of this set and its relations to the 
number of ancillary equilibria are discussed in [8], but we observe that in the case of the 2- 
cycles, C is of dimension 5. We have taken the square root of d for computational convenience. 
We denote by C the interior of C 

The rigidity matrix R of the framework is the Jacobian || restricted to the edges in E. We 
denote it by R = ff |#. The relevant definitions from rigidity theory are: 

1) Static rigidity: A framework is said to by statically rigid, or simply rigid, if for any d G C, 
there are only a finite number of frameworks, modulo rotation and translation of the plane, 
such that 6\e(p) = d. 

2) Infinitesimal rigidity: A framework is said to be infinitesimally rigid if there are no van- 
ishingly small motions of the vertices, modulo rotations and translations of the plane, that 
keep the edge-length constraints satisfied. This translates into llT6l : 

85 

rank( — \e) = 2n — 3. 
ox 

3) Minimal rigidity: A framework is said to be minimally rigid if none of the m frameworks 
with m — 1 edges obtained by removing one edge is rigid. 

B. Directed formation control 

We formalize in this section the type of control system considered in this work. We are given 
a graph G = (V,E) which is assumed to be minimally rigid with |V| = n and l^l = m. Let 
d G C be a feasible edge-lengths vector. The objective of the formation control problem is to 
find a decentralized control law, where the information flow is given by G, that will stabilize 
the system around a framework with inter-agent distances given by d. 

In more detail, each agent with position G M 2 is represented by a vertex Vi in V. The 
dynamics of Xi is allowed to depend only on the relative position of agents Xk for which there 
is an edge originating from and ending at t> fc : 

ii = Ui(x k , xl, . . .), where (vi, v k ), {v u vi), . . . G E. 
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In the case of directed control, it is easy to see that one cannot in general ask an agent to 
satisfy more than two edge lengths constraints ifTTl . For a subtler analysis of the situation, we 
refer to [fT8ll . From now on, we always assume that G has a maximum outvalence of two. 

Given orderings of the edges and vertices of G, we define 

Zi X k X\ 

if edge e, links nodes v k to v\. We define, with a slight abuse of notation, 

the error in edge length. 

We can thus write the set £ d for formation control problems as 

S d = {x e M. 2n \ei{x) = 0,for all edges in E}. 

An important feature of the formation control problem is that it is defined up to a rigid trans- 
formation of the plane: if x G M 2n is a framework of G, frameworks obtained by a rotation and 
translation of x — we write them as A ■ x, for A E SE(2), the special Euclidean group [{8J — 
are equivalent to x from a formation control point of view. As a consequence, we can assume 
without loss of generality that the agents measure only the relative positions of other agents. 

In most work on formation control, it is assumed that the agents' control law depends on the 
desired edge-lengths solely through the error in edge lengths e{. 

u = u(ei) = u(zfzi — di). 

We consider here a broader class of systems by allowing the Ui's to depend explicitly on the 
objective distances d. We distinguish this dependence of the control on d, which is a parameter 
by opposition to a dynamical variable, by using a semi-colon: u = u(di, zj Zi — di). 

We have: 

Xi tliidi) &i) Zi 

in case agent i follows a single agent and 

Xi = uu(d k , di] e k , e h z%zi)z k + u 2i (d k , d t ] e k , e h z£zi)zi 
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in case agent i follows two agents. For a complete justification of this model, see [fT9l . 

We established in [[HI a few conditions a feedback control law had to satisfy in order to yield 
a well-defined formation control system. We use them here to define a class of feedback control 
law for formation control: 

Definition 6. A feedback control law Ui is compatible with a formation control problem if 

1) Ui(dj] Cj) is such that Ui(dj\ 0) = if agent i has one co-leader. 

2) Ui(dj,dk;ej,ek,Zj ■ Zk) is such that ut(dj,dk]0,0, z) = for all z if agent i has two co- 
leaders. 

We accordingly define the class of controls U to be all twice differentiate control laws such 
that Ui(di] e») = and Uj(di, dj\ e*, ej, •) = have an isolated zero for e« = tj = 0, depending 
on whether the agent has one or two co-leaders. 

From the above discussion, we conclude that: 

Proposition 1. Let Ed be the set of design equilibria for a formation control problem with 
underlying information flow graph G. If the graph G is rigid, the set E^/SE(2) is finite. In 
other words, the set of design equilibria is finite up to rigid transformations. 




Fig. 7: The subgraphs with vertices {vi,v 2 ,v 3 } and edges {ei,e 3 ,e 4 } or vertices {vi,v 2 ,v 3 ,V4} 
and edges {ei, e 3 , e 4 , e 5 , e 7 } define subformations. The subgraph with vertices V = {v 1,1*4,^5} 
does not, for any set of edges, since outgoing edges point to vertices not in V. 



We end this section with the definition of subformation, illustrated in Figure |7} given a graph 
G = (V, E) underlying a formation, a graph G' = (V',E') underlies a subformation if G' is a 
subgraph of G where all outgoing edges from vertices in V are included in E '. 
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X 4 . 



Xi 



X-2 



•i'3 



(a) 



X i Xi/7 X 4 

/-L. x 2 

(b) 



^2 



•'V. 



(c) 



£4 

(d) 



x 3 



Fig. 8: Four formations in the plane that are not equivalent under rotations and translation and that have the same 
corresponding edge lengths, (a) is the mirror-symmetric of (c) and (b) is the mirror-symmetric of (d). 



C. The 2-cycles formation 



The 2-cycles is the formation represented in Figure lb Let Xi 6 



1 ... 4 represent the 



position of 4 agents in the plane. We define the vectors 



Z\ — X2 — X\\ Z2 — £3 — X2] Z3 — X\ — X3; Z4 — £3 — X4] Z§ — £4 — X\ 



(6) 



(7) 



Hence a general control law for such a system is 

±1 = u 11 (d 1 ,d 5 ]e 1 ,e 5 ,zlz 5 )z 1 +u 1 2(d 1 ,d 5 ;e 1 ,e 5 ,z'[z 5 )z 5 
x 2 = U2(d 2 ;e 2 )z2 
£3 = u 3 (d 3 ;e 3 )z 3 
X4 = u^d^e^jZi 

The set of design equilibria Ed for the 2-cycles is of cardinality 4, up to rigid transformations, 
since there are four frameworks in the plane for which e, = 0; they are depicted in Figure [8} 

In general, the set 8 a of ancillary equilibria depends on the choice of feedbacks Ui Eli. Due 
to the invariance and decentralized nature of the system, some configurations belong to S a for 
all elements of U: 

Proposition 2. The set £ contains, in addition to the equilibria in the frameworks charac- 
terized by 

1) Zi = for all i, which corresponds to having all the agents superposed. 

2) all Zi are aligned, which corresponds to having all agents on the same one -dimensional 
subspace in M. 2 . These frameworks form a three dimensional invariant subspace of the 
dynamics. 
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3) e 2 = e 3 = e 4 = 0, z\ and z 5 are aligned and so that 

ui(e 1 ,e 5 ,zlz 5 )\\z 1 \\ = ±u 5 {ei,er n zlz$)\\zr,\\, 
where the sign depends on whether Z\ and z 5 point in the same or opposite directions. 

This result is straightforward from an inspection of Equation ([7]). Frameworks of type 2 
above are not infinitesimally rigid and they define an invariant submanifold of the dynamics. 
Frameworks of type 3 appear whenever a vertex of the information flow graph has an outvalence 
of two. 



D. Singular formations for n = 4 agents. 

Let us gather the m vectors Zi G M 2 in z = {z\, z 2 , 



D(z) 







zl 
... 









... 
... 

zl 



o2m 



We set 



tmx2m 



Hence D(z) G 

In order to use Sotomayor's theorem, we single out a particular type of frameworks which, 
even though they are infinitesimally rigid, show a certain degree of degeneracy. We first observe 
that, in general, the angle between z\ and z 5 is not uniquely determined by the edge lengths. We 
define S to be set of edge lengths such that at least one of the four frameworks corresponding 
to d has zi parallel to z 5 with the notation of Figure [9] 

S = {d E £ s.t. zi parallel z 5 for one framework at least.} 



and 5 = 5n C . 

We will need the following properties of S: 

Lemma 1. The following properties of S hold: 
1) S is of codimension one is C 
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X 2 




Fig. 9: The formation in (a) is such that (||«i||, ■ ■ • , H^sll) ^ S, whereas (||zi||, ■ • • , \\z5W) € S for the formation 
depicted in (b) 

2) The formations corresponding to edge lengths in So are infinite simally rigid. 

Proof: For the first part, observe that we can parametrize S by first choosing a feasible 
di,d 2 ,d 3 yielding a triangle xi,x 2 ,x 3 and one additional parameter giving the signed length of 
z 5 , with the sign referring to z 5 going in the same direction as z\ or the opposite direction. 
Because the set of feasible d\,di,d$ is included in £, we see that we need 4 parameters to 
describe a formation in S and hence it is of codimension one. 

For the second part, we have that the rigidity matrix of the 2-cycles is given by R: 



T 

4 


T 

-4 











4 


T 

-4 





T 

-4 





T 

4 











-4 


T 

4 


T 

4 








-4 



Some simple algebra shows that one has 

R = D(z)A$ (9) 

where we recall that A m is the mixed adjacency matrix. In the case of the 2-cycles, the mixed 
adjacency matrix A m G IR 5x4 is of rank 3. The cokerneQ of A m is spanned by [0, 0, 1, 1, 1] T 
and [1, 1, 1, 0, 0] T . Hence, the cokernel of A$ is four dimensional and spanned by the vectors 

2 The cokemel of a linear map / : A — ¥ B is the quotient space _B/im(/). Its coimage is A/ker(f). 
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[0, 0, 1, 1, 1] T <g> [1, 0] T , [0, 0, 1, 1, 1] T <g> [0, [1, 1, 1, 0, 0] T ® [1, 0] T and [1, 1, 1, 0, 0] T ® [0, 1] T 
The matrix -D(z) is of full rank unless Zi — for some i, which corresponds to two agents 
superposed. We thus have that D(z) is of full rank for formations in So. Because D(z) is of 

(2) 

full row rank, R is of full (row) rank if A m ' projects onto the coimage of D(z), which is easy 
to see by inspection. ■ 

VI. The 2-cycles is not almost surely stabilizable 

A. Statement of the main result 

The main result of this paper is to show that given any robust control law, there is an open 
set of target configurations that are not almost surely stabilizable for the 2-cycles; equivalently, 
to show that for any control law that robustly stabilizes a design equilibrium, there is an open 
set of d's in C for which an ancillary equilibrium is stable. 

Theorem 2. The 2-cycles formation is not robustly almost-sure stabilizable for an open set of 
design frameworks. 

The result readily extends to formations that contain the 2-cycles as a subformation: 

Corollary 1. Any formation that contains the 2-cycles as a subformation is not robustly almost- 
sure stabilizable for an open set of design frameworks. 

B. Proof of the main result 

We rewrite the dynamics in terms of the z variables 



Zl 


= U 2 Z 2 - UiZi 


- u 5 z 5 


z 2 


= u 3 z 3 - u 2 z 2 




z 3 


= UiZi + u 5 z 5 


- u 3 z 3 


ZA 


= U 3 Z 3 - U4Z4 




Zb 


= U4Z4 — U\Z\ 


- u 5 z 5 
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where the dependence of the u { on di and Zj's is given in Equation ([7]). We denote by F(z) the 



right-hand side of Equation ( fT0| ) and set 

7" = | Ui G W}. 



We now show that the system of Equation 10 has the logistic equation as normal form at S 



Theorem 3. A transcritical bifurcation at frameworks with d in an open set in Sq is generic for 
systems in T. 

We will prove Theorem [2] as a corollary of this result. We prove Theorem [3] in several steps. 



Proposition 3. Let d G Sq. There is a non-zero vector w G IR 10 such that w T ff L=o,d 

dFi 
dd 



w T ^K\d = for at least one framework attached to d. 



We denote by u x (di, <i 2 ; x, y, z) the derivative of u with respect to x, and similar definitions 
hold for u y and u z . We need the following lemmas: 

Lemma 2. Set 

z\ = 2(u lx Zi + u 2x Zj) 

and 

z'j = 2(u ly Zi + u 2y Zj) 
if Zi originates from vertex with outvalence two and outgoing edges Zi and Zj, and 

Z^ 2xi x Zi 

if Zi originates from a vertex with outvalence one. The Jacobian at a design equilibrium of a 
formation control system is 

d -f = A^D{z'YD{z). (11) 
oz 
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Proof: We first observe that 





.^ux{ei,ej,zjzj)zi = 2u lx Zizf + ul 2 + 2u lz Zizj 



dzj 



ux(ei,ej,Zi ZjjZi = 2u ly z i z j + 2u\ z z i z j 



■^-u 2 (e ij e j ,zfz j )z j = 2u 2x z j zf + 2u 2z z j zf 
-^-u 2 (e i ,e j ,zjz j )zj = 2u 2y ZjzJ + 2u 2z zjzj 



3 



where we omitted the arguments of the functions on the right-hand side. Recall from Definition [6] 
that u z = u = at a design equilibrium. Hence, if Z\ originates from a vertex with two outgoing 
edges with 

Zi = Fi(z) = . . . - ui(e h ej, zjzj)zi - u 2 (e u ej, zjzj)zj 



then: 



dz. 



-2u lx ZizJ - 2u 2x Zjzf = -z[zj 



Similarly, 



^ = -2z'z T 
dz, i > ' 



"J 

If z k originates from an agent with a single leader, we have: 

9Fk t IT 

— - u xZj z 3 - z 3 z 3 . 

In general, if z k appears in F h then 

dz k - ±ZkZk 

where the sign is negative if both z\ and z k are leaving the same vertex and positive if z h is leaving 
the vertex to which z k leaves. Putting the equations above together and recalling Definition [5] of 
A e , we get the result. ■ 
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We define 



„'/ _ J \ dm „ , du 5 
z l - z l + ddi Zl ddi Z5 



, du2 „ 
Z 2 + ddk Z2 

4 + If ^3 
Z 5 "T ad 5 1 dd 5 Z 5 



< 4' 



(12) 



Lemma 3. The Jacobian of F with respect to the parameters d at a design equilibrium is given 
by 



dF 
~dd 



A®D(z") T . 



(13) 



Proof: We have 



dF 1 dm du 5 

dd- 1 =~ i dd; + Ulx)zi ~ i dd- + U5x)z5 



and similar relations for the other entries Some algebraic manipulations yield the result. 



Lemma 4. Let d E S , then w is a left eigenvector of ^\ ei =o,d with eigenvalue if and only 
ifw Td ^U =0 , d = 0. 

Proof: We claim that when z\ is parallel to Z5, we can find invertible diagonal matrices D\ 
and D 2 such that 

D(z') = D x D{z) = D 2 D{z"). (14) 

Indeed, for i = 2,3,4, it is immediate from the definitions of Zi, z\ and z" that there exists 
afii 7^ such that 

Zi = aiiz'i = (3iz", i = 2,3, 4. 

The (Xi,/3i are the entries of D 1 and D 2 . From Lemma [2} we have 

dF 



dz 



A^D(z') T D(z) 



(15) 



Putting Equations (13), (14), (15) together, we obtain that 



^F nn(^ dF 
-^D 2 D(z) = —. 
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If all the Zi are non-zero, then D 2 D(z) is of full rank and we conclude that 

w T — = & w T — = 0. (16) 

od az 



Formation control systems are invariant under an action of the Euclidean group SE(2) on 
M 2 . Hence the Jacobian of the dynamics in the x variables will always have at an equilibrium 
three zero eigenvalues corresponding to the three dimensions of SE{2). For the dynamics in 
the z variables, there are additional zeros from the redundancy of the z. For example, described 
in the z coordinates, the 2-cycles has 10 dimensions, but four degrees of freedom are lost to 
^1 + ^2 + ^3 = z 3 + z^ + z 5 = and one additional degree of freedom is lost to the SE{2) 
invariance, since the invariance under translation is taken into account in the z variables and 
only the invariance under rotations, which is one dimensional, remains. The following result 
addresses this point. 

Corollary 2. Let G be the graph of a minimally rigid formation with edge adjacency matrix A e . 
The eigenvalues of the Jacobian ofF(z) at a non-singular design equilibrium are the eigenvalue 
zero with algebraic multiplicity 2n — 3 and the eigenvalues of 

J = D(z)A e D(z') T . (17) 

Proof: The result is a consequence of Theorem 1.3.20 in [|20l applied to Lemma |2j ■ 
In the remainder of the paper, whenever we refer to eigenvalues and eigenvectors of F, we 
will refer to the eigenvalues and eigenvectors that do not correspond to the redundant description 
of the system, and to which we have access thanks to Corollary [2} 

Corollary 3 (Singular formations). Let d E S . The Jacobian of the 2-cycles formation is 
generically of corank 1 for at least one framework attached to d. 

Proof: A direct computation using Corollary [2] and the edge- adjacency matrix of the 2-cycles 
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gives 



J = D{z)A e D{z 



l\T 



-z x Z l 


z l z 2 


U 


n 
u 


z 1 z 5 





T I 

~ z 2 Z 2 


T l 

z 2 z 3 








T l 

Z 3 Z-y 





T l 

Z 3 Z 3 





T l 

Z 3 Z 5 








T i 

z 3 z A 


T i 

z^ Z^ 





T I 

z 1 z 5 








T i 

z A z 5 


T i 

z b Z 5 



By Corollary [4] in the appendix, u[ ^ generically when it vanishes, which correspond to 
design equilibria. Hence z[ are generically non-zero. For the framework attached to d E So such 
that Zi is parallel to z 5 , the first and last column are multiples of each other, and it is easy to 
see that the first four columns are linearly independent. The corank is higher if, in addition, one 
of the Zi is zero. ■ 

We now prove Proposition [3] 

Proof of Proposition ^ Consider a framework with d E S and z x parallel to z 5 . From 
Corollary [3J we know that ^ is generically of rank 4. Let w be an eigenvector corresponding 
to the zero eigenvalue. We conclude using Lemma [4] that w T ^ = 0. ■ 

For clarity of the exposition, we now restrict ourselves to the system 



±i = u(d 1 ;e 1 )z 1 +u(d 5 ;e 5 )z 5 
x 2 = u(d 2 ;e 2 )z2 
x 3 = u(d 3 ;e 3 )z 3 
= u(d 4: ; e^z^ 



(18) 



The preliminary results, proved in greater generality, make the extension of the proof below to 
the more general system easy. 

Proof of Theorem^ Fix do = (di, d 2 , d$, d±, d^) E Sq. We consider the one parameter 
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system where only juelis allowed to vary: 

±i = u(dx;e 1 )z 1 + u(d 5 ;e 5 )z5 
x 2 = u(d 2 ;e 2 )z 2 

(19) 

x 3 = u(d 3 + fi;zjz 3 - (d a + ^))z 3 
X4 = u(d^e4)z4 

and the corresponding equations in z variables. We prove that conditions (1), (2) and (3) of 
Theorem [T] are generic for F. From Corollary [3] and the fact that vl 7^ generically at a zero of 
u (Corollary [4]), we know that the Jacobian of the 2-cycles at S has a unique zero eigenvalue 
(that is not a result of the redundant description of the system) generically for F G F. Hence 
condition (1) is verified. Condition (2) follows from Proposition [3] 

Condition (3) takes the following form: for w and v left and right eigenvectors of ^ 
respectively, 



T d 2 F \ (9 2 Fj 



where w and t> depend on the design equilibrium. Using the relations established in Lemma [2] 
for partial derivatives of Fi, the triple sum in the above Equation can be explicitly evaluated. 
Rearranging terms and introducing the constants c\, c 2 (which depend on the design equilibrium), 
we obtain a linear combination of first and second derivatives of u: 

d 2 F 

wT ^d^^ v ' ^ = ° lUx + C2Uxx ' 

In order to verify that there is an open set of design equilibria for which the functions C\ and 
C2 are non-zero, since the functions are continuous, it suffices to find one design equilibrium 
at which it is the case. It is easily verified, for example, at the design equilibrium with edge 
lengths corresponding to the framework with x\ = (0,0), x 2 = (—2, l),x 3 = (0, —1) and x 4 = 
(1/2,-1/4). 

On that open set, we thus have w T ^-(v, v) — only if 

c\u x + c 2 u xx = 

when u vanishes (since, by definition, u vanishes at design equilibria). Let C C J 2 (IR, M) be 
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defined by the equations 

c\u x + c 2 u xx = 

and 

u = 0. 

C is of codimension 2 in J 2 (IR, E) whereas the image of the 2-jet extension of w is of dimension 
1. Hence it is transversal to C if and only if u ^ or C\U X + c 2 u xx ^ 0. We conclude using 
Thorn's transversality Theorem [4] that w T ^-(v,v) is generically non-zero. 

Using a similar reasoning as above, we can conclude that w T J^t> is generically non-zero. 

■ 

Proof of Theorem |2| - 

We illustrate the situation in Figure [TOj We will show that there is a set of positive measure in 
C which cannot be made robustly almost surely globally stable. We do so by showing that for 
any framework attached to distances in that set, the existence of a stable ancillary equilibrium 
is generic for T . 

Denote by S £ a tubular neighborhood of S: 

S £ = {deC s.t. 3 d e S with \\d - d \\ < e} 

and S £ = S £ nSo- The set S £ contains frameworks where z\ and z 5 are close to parallel. These 
frameworks are infinitesimally rigid and non-singular. Let d E S £ and d e S be such that there 
is — e < [i < e with d = do + (0, 0, /x, 0, 0). Such do and \i exist by definition of S £ . 

Because the system is invariant under mirror symmetry J8), the stability properties of the 
equilibria (a) and (c), (b) and (d) in Figure [8] are the same. Assume without loss of generality 
that u is such that the design equilibria for the frameworks with x 2 and £4 on different side of zj, 
are stable. Because for a robust u, the system undergoes a transcritical bifurcation when ji = by 
Theorem [3j and because u' 7^ generically when u vanishes, we have that for e small enough, 
8 a contains the framework where z\ is parallel to z 5 for all frameworks with — e < /x < e. 
Furthermore, for either /x > or /x < 0, we have that this framework is asymptotically stable, 
i.e. S s H £ a 7^ 0. Hence, there is a set of positive measure of target frameworks in Sq which 
contains a stable ancillary equilibrium and thus the system is not robustly almost-sure globally 
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x ids 
^2 di+ii-« 




Fig. 10: We illustrate the stability properties of ancillary and design equilibria around Sq. Let the vector 
(di,d2,d3,d^,ds) € Sq. The horizontal dashed line corresponds ancillary equilibria and the slanted line that 
intersects it to design equilibria. They coincide at fi = 0, as seen in Proposition V-C for /i f^ configurations in So 
are ancillary equilibria. For \i\ < 0, there is an ancillary equilibrium with 62, 63, = but e\ = £\ and e$ = £2 
and Z\ and z 5 aligned. It is illustrated in the top-left corner of the figure. This equilibrium is moreover stable. For 
H2 > 0, there is a similar ancillary equilibria with e\ — £3 and es = £4, illustrated in the bottom-right corner, but 
this equilibrium is unstable. We see that around the bifurcation value fj,Q, there is a transfer of stability from £ 4 to 
£ a . The orientation may be reversed (i.e. /ii > 0, H2 < and all else the same in the figure) depending on the sign 
of the second derivatives in Theorem [3] 



stable. 
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Appendix 

The main tool handling genericity and robustness in function spaces is Thorn's transversality 
theorem. We will arrive at the result by building onto the simpler concept of transversality of 
linear subspaces. 

Let A,B dW 1 be linear subspaces. They are transversal if 

R n = A © B, 
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Fig. 1 1: We let V be the property of vanishing with a zero derivative. Then 1 — V is generic and thus V is not 
robust. Let u(x) be a function which satisfy V. For almost all perturbations, it will either vanish with a non-zero 
derivative — as illustrated with u{x), dashed curve — or not vanish at all — as illustrated with u(x), dotted curve. 
Both u(x) and u(x) are transversal to the manifold defined by y = everywhere, whereas u(x) is not. 



where © denotes the direct sum. For example, a plane and a line not contained in the plane are 
transversal in IR 3 . The notion of transversality can be extended to maps as follows: given 



F 1 : W 1 ->■ R m and F 2 : R l -> 



5 ;/( 



we say that Fi and F 2 are transversal at a point (xi, £2) G E n x IR' if one of the two following 
conditions is met: 

1) F 1 (x l ) ^F 2 {x 2 ) 



2) If Fi(xi) = F 2 (x 2 ), then the matrix 



dFx dF 2 
dx dx 



1 T 



is of full rank. 



In particular, if I + n < m then F% and F 2 are transversal only where they do not map to the 
same point. 

The notion of transversality that is of interest to us is a straightforward extension of the 
transversality of maps: 

Definition 7 (Transversality). Let F : M N be a smooth map and let C be a submanifold of 
N. Then F is transversal to C at a given point if, at that point, F is transversal to the embedding 
i:C^N ofC into N. 

The definition is best understood from the following example: 

Example 3. Take N = IR 3 with coordinates u,v,w and C be the u-v plane. Let F : IR — >■ IR 3 : 
x — > [x, 2x, 3x] T . Then the map F is transversal to C everywhere since F(x) £ C for i ^ 
and at x = 0, ^ is not in the tangent space of C. 
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Let F, G : M — > N be smooth maps between smooth manifolds M and N . We say that F 
and G are O-equivalent at x if F(x ) = G(xq), 1 -equivalent if in addition to being O-equivalent, 
ffUo = f^Uoi an d so forth. We define the k-jet of a smooth map to be its k-equivalence class: 

Definition 8. The k-jet of F : M N at x is 

J x ( F ) = {G : M ^ N s.t. G is k-equivalent to F}. 

Hence, the 0-jet of F at x is F(x ); the 1-jet is (F(x ), ^\ Xo ), etc. For example, the constant 
function and sin(x) have the same 0-jet at x = and x and sin(x) have the same 1-jet at 0. 
We define: 

J k {M, N) = Space of k- jets from M to N. 

Given a function F : M — > N, we call its k-jet extension the map given by jp(x) : M — > 
J K (M,N):x^(F(x),^(x),...,^(x)). 

Example 4. Let M = N = R. The k-jet space is J k (R, R)=RxRx...xR = R k+2 . Take 
F(x) = sin(x); the 2— jet extension of F is 

Jsin( x ) = (z,sm(x),cos(x),-sin(a;)). 

If we take M = N = R 2 and F(x) = Ax for A e R 2x2 , then j k Ax (x) = (x, Ax, A, 0, ... , 0). 

While to any function F : M — >■ N , we can assign a k-jet extension jp : M — > J k (M, N), the 
inverse is not true: there are maps G : M — >■ J k (M,N) which do not correspond to functions 
from M to iV as there are some obvious integrability conditions that need to be satisfied. For 
example, if we let 

G : R n -)• J^R") : = (x, Ax, 5), 

then G is a 1-jet extension of a function if and only if B = A. 

The power of the transversality theorem of Thorn is that it allows one to draw conclusions about 
transversality properties in general, and genericity in particular, by solely looking at perturbations 
in jet spaces that are jet extensions — a much smaller set than all perturbations in jet-spaces. 

We recall that the C r topology is a metric topology. It is induced by a metric that takes into 
account the function and its first r derivatives. We have: 
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Theorem 4 (Thorn's transversality). Let C be a regular submanifold of the jet space J k (M, N). 
Then the set of maps f : M — >■ N whose k-jet extensions are transversal to C is an everywhere 
dense intersection of open sets in the space of smooth maps for the C r topology, 1 < r < oo. 

A typical application of the theorem is to prove that vector fields with degenerate zeros are 
not generic. We here prove a version of this result that is of interest to us. 

Corollary 4. Functions in C°°(M) whose derivative at a zero vanish are not generic. 

In other words, the corollary deals with the intuitive fact that if u(x) = 0, then generically 

u'{x) ^ 0. 

Proof: Consider the space of 0-jets J°(M, M). In this space, let C be the set of 0-jets which 
vanish, i.e. C = (x, 0) C J°. A function u is transversal to this set if either it does not vanish, 
or where it vanishes we have that the matrix 

1 1 

^ 

dx. 

is of full rank. Hence, transversality to C at a zero implies that the derivative of the function is 
non-zero. The result is thus a consequence of Theorem |4j ■ 
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